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FOREWORD 


The work described herein was performed at Mechanical Technology Incorporated 
(MTl) under NASA Contract NAS3-15334 with Mr. Robert E. Cunningham, NASA- 
Lewis Research Center, as Project Manager. At MTI, Dr, Robert H. Badgley, 
Manager, Machinery Dynamics Center, was the Program Manager. The work 
reported is part of a continuing program whose objectives are, firstly, to 
develop, through analytical and experimental investigations, elastomer 
dynamics technology and, secondly, to reduce this technology to practical 
designer-oriented form, so that it may be used to calculate the dynamic 
mechanical properties of elastomeric damping materials. 
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I SUMMARY 


This report documents an experimental investigation to determine the dynamic 
properties of elastomer test sauries at constant temperature. The objective 
of the investigation was to identify the influence of excitation frequency, 
preload, shape, and deformation amplitude on the dynamic properties under shear 
and compressive loading. 

An existing forced-vibration resonant -mass-type test rig was refined to con- 
trol the ambient temperature around the elastomer san5>le and monitor the power 
dissipated in the specimen. These two constraints ensure a constant thermal 
state of the elastomer. Polybutadiene, a broad temperature range (BTR) 
elastomer, was selected as the test material for this investigation. 

Three compression specimens with varying shape factors were tested at preloads 

corresponding to a compression of 2-1/2 and 5 percent of the initial specimen 

thickness. The ambient temperature and power dissipation were monitored and 

controlled effectively such that the elastomer temperature was 32®C to within 
0 

1 C. The test frequencies were varied from 100 to 1000 Hz. Test results are 
presented as complex stiffness values for each specimen as a function of fre- 
quency. Changes in stiffness as a function of frequency or preload were found 
to bt small for the compression specimens, and within the test frequency range 
a simplfc pwer law relationship was used to represent the variation with fre- 
quency. 

Three shear specimens were tested at preloads corresponding to nominal i. trains 
of 0, 2-1/2 and 5 percent. The frequency band was also 100 to 1000 Hz. The 
complex stiffness in this case was somewhat higher than that observed with the 
compression specimens, but the difference is consistent with the difference in 
static stiffness as calculatea by conventional shape factor formulas. Varia- 
tions of complex stiffness with either frequency or preload were found to be 
small, and a sinq>le power law relationship was again used to represent the 
variation with frequency. 

The compression samples, although covering a four to one range of dlameter-to- 
height ratio, were each designed to have similar stiffnesses by varying the 



number of Individual elastomer elements, and applying published shape factor 
values. The measured dynamic stiffnesses of the three compression samples were 
similar. 

The shear samples were each designed to have similar stiffnesses by varying 
the number of elements in Inverse proportion to the sheared area. The resultant 
dynamic stiffnesses were similar for each of the shear samples. 
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II. INTRODUCTION 


Elastomeric dampers have been successfully applied in the control <> uk,i. 'leal 
vibration problems. Increasingly demanding applications, particu , arly to ibra 
tion problems of rotors and power transmission shafting, arc being considered. 
However, the availability of design-oriented data for such application': Ls 
limited. Dynamic testing under accurately controlled conditions is the means 
to develop this data and to determine the influence of in^ortant geometrical, 
environmental, and chemical design parameters. Physical properties of 
elastomeric materials can then be evaluated and used in the selection of 
materials for particular applications. 

Substantial advances have been reported in the development of experimental 
methods for dynamic testing, and the determination of the influence of chemical 
composition on the observed dynamic response of the material. Chiang, et. al., 
in Part I of this series [1]^, described procedures for testing elastomer 
specimens in a forced-vibration, resonant mass test rig. The ability of the 
resonant mass approach to obtain data over a wide range of frequencies, 
amplitudes, and preloads was demonstrated. In addition, the significance of 
uncertainties in temperature variation was identified. 

A narrative account of various test methods and some of their refinements 
developed during the past two decades has been recently presented by Miller 
[2]. In addition to the measurements of frequency and amplitude of deforma- 
tion, a major concern is with the temperature level and gradient in the elasto- 
mer. 't has been suggested that the thermal state of the material is one of 
the prime candidates dominating its dynamic behavior. Hence, monitoring the 
temperature of operation and the energy dissipated in the elastomer becomes 
essential in an experimental setup for elastomer testing. 

The significance of dynamic testing in comparison with simple static tests has 
been well established. Larsen [3] has demonstrated that the dynamic stiffness 


Numbers in brackets <\>note References listed at the end of this report. 
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is generally higher than the static spring rate. The temperature dependence 
of the two properties as reported by Larsen is even more Interesting In that 
although the dynamic stiffness decreases with an Increase In cemperature, the 
changes In the long chain molecules as a result of higher temperature results 
In Increased static stiffness. 

''.ardlllo [ 4 ! used a force vibration type hydraulic tester ~o determine the 
stiffness and damping properties of a number of materials. The effects of 
frequency, amplitude, preload, and temperature were also considered. In ;:he 
case of butyl stock. It was shown that Increasing the carbon black content In 
the compound results In Increased dynamic stiffness and damping characteristics. 
Heyer and Sommer [Sl have considered con^oimds based on styrenebutadlene and 
polybutadiene elastomers In their Investigation with regard to the effects of 
compositional and testing variations. Their results concerning the Influence 
of carbon black content are In agreement with those reported by Cardlllo. In 
a more recent work, Sommer and Meyer [6] have summarized the effects of varia- 
tions In rubbers, fillers, oils, level of crossllhklng and processing ^ "ors 
such as mixing, curing, and storage on the dynamic properties of e' ^rlc 
products. 

The work reported herein was motivated by a need to further liq>rcve the dynamic 
test methods for elastomers with particular reference to thermal conditions. 

By controlling the ambient teBq>erature and the power dissipated In the elastomer 
specimen, the uncertainties In the dynamic behavior of the test material are 
reduced to a minimum. The primary objective of the Investigation Is to Identify 
t) e Influence of excltatlc.i frequency, preload, shape factor, and deformation 
amplitude on the dynamic properties of elastomer specimens under conditions 
wherein the thermal state of the material remains unchanged. 

A previously designed [1] forced vibration type of nesting apparatus was 
modified to control the tenq>erature of the two metallic plates holding the 
elastomer san^le. Also, heating and cooling systems were designed to control 
the ambient temperature In a closed chaoiber In «dilch the test specimens are 
located. An electronic circuit to monitor the power dissipation was designed 


with the objective of limiting the temperature gradients in the elastomer 
specimens. Polybutadiene, being a known broad temperature range elastomer, was 
selected as the test material. The influence of excitation frequency, amplitude 
of deformation, preload, and shape factor were examined um'er a constant ambient 
temperature and very low and bounded temperature gradients In the test specimen. 
Both shear and compression loadings were considered with suitably designed 
specimens. The compression speclm* i were designed with varying shape factors 
with the objective of identifying the variations in the dynamic behavior of 
the elastomer as a function of specimen geometry under isothermal conditions. 

The experimental details of the apparatus are described in detail in the next 
section, which is followed by a fairly comprehensive analysis of the dynamic 
response of the test rig. Experimental results for both compression and shear 
samples are then presented. A coiq>lete summary of these results is presented 
in a tabulated form. Finally, the dynamic behavior of the elastomer . s a 
function of the experimental conditions described above is discussed. 



Ill, EXPERIMENT DESIGN AND TEST PROCEDURE 


A forced-vibration resonant mass type of apparatus was used to test the dynamic 
properties of the elastomer samples, A general description of the test rig is 
presented below before discussing the design characteristics. The details of 
selected material, test specimens, required instrumentation and the test pro- 
cedure are also presented in this section, 

A, Description of Elastomer Test Rig 

The basic features of this test rig have largely been developed in the course 
of work conducted under a previous phase of the present NASA Contract. The 
test rig was designed to Impose precisely measured uniaxial vibration ampli- 
tudes upon a selected elastomer sample at desired test frequencies and under 
selt:ted force preloads. 

The rig is capable of testing, through the base-excitation resonant-mass 
technique, elastomer samples of virtually any size and shape (test sample 
cavity is a cylinder approximately 12.7 cm (5 in,) high by 12,7 cm (5 in,) in 
diameter over a wide frequency range. Selection of the resonant mass to 
match elastomer sample properties permits testing at very high amplitudes at 
resonance, with correspondingly lower amplitudes at off-resonance conditions. 

Thus, maximum test amplitudes are limited by elastomer characteristics and 
shaker input power rather than by the test rig itself. 

The vibration input to the test rig in all tests was obtained from a commercial ly- 
available electromagnetic shaker system^' capable of delivering 66,700 N (15,000 
lb) for :e in the sinusoidal mode of vibration. 

Prior to the acquisition of the elastomer test data reported here, ambient 
temperature control facilities in the form of heating elements and cooling 
coils embedded in the sample base holder and circuitry to monitor the power 
dissipated in the specimen were added. Uncertainties due to temperature 


Ling Electronics, Model 335A Shaker With PP-35/70 VC Power Amplifier and 
SCO-100 Servo Control Center. 



variations are now minimized by accurately controlling the ambient tempera- 
ture, and by limiting the heat generated in the specimen. In addition, new 
holding plates for the elastomer samples were made; two different configura- 
aons provided shear and compression loading of the samples, respectively. 

B. Design Characteristics of the Test Rig 

A schematic of the elastomer test rig is shown in Figure 1 and a layout 
drawing of the same rig in Figure 2. Figure 3 shows the complete rig 
mounted on the shake table. 

The test rig was designed to meet the following functional requirements: 

1. Very low residual damping. 

2. An inertia mass loading of the elastomer test sample with weights 
ranging from approximately 0.8 to 227 kg (1.8 to 500 lb). 

3. A means for accurately controlling force preloading of the elastomer 
test sample. 

4. Temperature control of the ambient air within the test specimen 
cavity. 

The rig employs a unidirectional vibrating mass-spring system. In its 
simplest form, this concept might ^.onsist of an elastomer spring to which 
various amounts of inertia mass are bonded and the system could be completely 
free of residual damping or stiffness contributions from test rig elements. 
However, the practical requirements of balancing rather large inertia masses 
(up to 227 kg. (500 lb))on top of the elastomer spring and the desired mechanism 
to provide elastomer sample preloading at all test frequencies introduce addi- 
tional elements which have finite stiffnesses and small amounts of parasitic 
damping. 

Specifically, two radial guide bearings for the inertia mass and two air 
cylinders for the preloading of the elastomer samples are required. The 
stiffness and damping properties of these elements have an influence on the 
test rig response, but as described in Section IV, these properties are indepen- 
dently measured and accounted for. 
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The resonant mass bonded to the top of the elastomer specimens consists of 
a number of elements. For high frequency tests, in which only a minimum of 
mass is required, the resonant mass can be reduced to the top plate of the 
specimen holder and the preload piston. The top plate is made out of aluminum 
and the preload piston is made out of titanium for minimum weight. Increased 
mass for medium or low frequency testing is provided by rigidly attaching a 
long rod to the top of the preload piston, and adding steel weights. The 
weights are centered by the rod and are axially restrained by spacers of 
various lengths and a locknut near the upper end of the rod. 

The rod itself receives radial support from two frictionless guide bearings 
(Figure 4 shows one of the dismounted bearings), each of which consists of a 
hub and 12 steel spokes. The hub is fitted over the end of the rod and is 
axially clamped to it. This kind of bearing arrangement provides good radial 
stiffness for reasonably high spoke tension, but provides only limited freedom 
for axial motion if overstressing of the spokes is to be avoided. Consequently, 
changes in preload on the test specimen must be accompanied by lowering or 
raising of the outer frame (to which the guide bearings are attached). 

Preload of the elastomer test specimens is controlled by a pair of air cylinders. 
The lower cylinder is located directly above the test specimen. Its piston is 
attached to the specimen's upper holding-plate and exerts a downward force 
on the specimen. The upper cylinder is located at the top of the test rig and 
its piston exerts an upward force on the resonant mass. When the resonant 
mass is large, the desired preload is usually achieved by pressurizing the 
upper cylinder and partially supporting the weight of the mass. When the 
resonant mass is small, the desired preload is usually achieved by pressurizing 
the lower cylinder and exerting a downward force in addition to the weight of 
the mass. 

The upper and lower preload pistons are sealed in their cylinders by two 
roiling diaphragm seals (Bellofram 3C-600-37-FPJ and 3-119-119-CBJ) . The air 
cylinder inlet hole is 0.762 mm (0.030 in) in diameter, making each cylinder 
essentially a "closed" cavity under vibration conditions. A constant preload 
force on the test specimen is maintained through regulation of air pressure in 
the cylinders. 
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The ambient temperature within the elastomer test specimen cavity is controlled 
through the combined effect of electric heating cables and water-cooled tubes 
installed in the base plate below the test specimens (see Figure 5 for a 
photograph of the installed heating-cooling elements) . The grooves in which 
the heating cable and the tubing are embedded have a reverse taper to ensure 
secure and rattle-free operation in a vibration environment. The electric 
heater is automatically controlled by a temperature controller, which obtains 
the necessary feedback from a thermocouple embedded in the metal of the base 
plate. 

Details of the temperature control system and of the ter..perature recording 
system are shown in Figure 6. 

Cj; rest Material Selection 

As discussed above, the temperature variations are minimized by accurately 
controlling the ambient temperature and the heat generated in the specimen. 

To reduce the effects of any temperature variations which remain, a Broad 
Temperature Range (BTR) elastomer was selected as the test material. BTR 
elastomers are claimed to show reduced variation of material properties with 
temperature. 

Two of the BTR-type elastomers under consideration for testing were silicones 
and polybutadiene. The final choice was polybutadiene, because of the 
difficulty, according to the manufacturer, of producing batches of silicone 
materials with consistent material properties. As with most elastomers, 
polybutadiene may be produced with various material hardness and damping rates, 
according to proprietary manufacturers' formulas. The polybutadiene chosen 
for the test described in this report carries the manufacturer's (Nichols 
Engineering, Inc., Shelton, Connecticut) designation NE X15SG, with a nominal 
hardness of 70 durometer (or Shore A Hardness) and 21 percent of "critical 
damping". This material had the highest hardness and damping of those avail- 
able from the manufacturer. 

The individual polybutadiene test specimens were bonded with a cyanoacrylate 
adhesive (Eastman 910 HMP) to their respective aluminum support plates. The 
choice of bonding agent was based upon the manufacturer's recommendations and upon 
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sample tests, which showed the epoxy (HYSOL AS7-4323) which was previously used 
for the bonding of neoprene and urethane elastomer specimens to steel and 
titanium surfaces yielding unsatisfactory results. It was found that successful 
bonding with Eastman 910 HMP is dependent upon smooth bonding surfaces, a thin 
film of adhesive and rapid application of a considerable clamping force to the 
bonded surfaces, 

D. Shear Test Samples 

The shear test specimens were molded as sheets 15 cm by 15 cm (6 in. by 6 in.) 
with a thickness of 3.2 mm (0.125 in. ), Individual test specimens of the following 
quantities and sizes were knife-cut from these sheets: 

1. Four strips - 25.4 mm high x 48.8 mm long (1.0 x 1.92 in.) 

2. Eight strips - 12.7 mm high x 48.8 mm long (0.5 x 1.92 in.) 

3. Twenty strips - I.l mm high x 48.8 mm long (0.2 x 1.92 in.) 

In the case of the first sample configuration! the four individual strips 
were glued with one of their large faces to the sides of a square aluminim 
block and with their other large face each to one additional aluminim holding 
block (see Figure 7 for a photograph of the assembly of the four-strip shear 
test sample). In the assembly of the test fixture, the four holding blocks 
were bolted to the shake-table, while the elastomer-supported center block 
was attached to, and tnerefore became part of, the resonant mass. 

The second shear test sample in its assembled configuration (Figure 8) was very 
similar to the first one, except that each of the assemblies was only of half 
the height. The two assemblies, placed on top of each other, then fitted into 
exactly the same space in the test fixture as the first test sample. The 
eight strips, deflecting in unison, were designed to have the same stiffness 
against deflection of the center mass relative to the mounting blocks, as the 
four strips in the first assembly (see Appendix A). 

The convenience associated with the assembly of individual four-specimen sub- 
assemblies becomes quite apparent upon inspection of the third shear test 
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Test Assembly of Four Elastomer Shear Specimens, Eacl 




sample (Figure 9). To maintain the same stiffness of the assembly with shear 
strips only 5.1 mm (0.2 in.) high, a total of twenty were needed. Since the 
individual strips were required to be free of interference by the strips above 
or below it, they were glued to aluminum blocks 6.4 mm (0.25 in.) high. When 
assembled in a stock of five, there is a height separation of 1.3 mm (0.05 In.) 
between the individual strips. 

E. Compression Test Samples 

Individual compression test specimens were shaped as c; llnders of various 
heights, rhe cylinders were molded to a uniform height of 12.7 mm (0.5 in.) 
and later cut and ground to the desired height. From individual specimens 
the following three test samples were assembled: 

1. Thirty cylinders - 12.7 mm diameter by 12.7 mm height (0.5 x 0.5 in.) 

2. Ten cylinders - 12.7 mm diameter by 6.4 mm high (0.5 x 0.25 in.) 

3. Three cylinders - 12.7 mm diameter by 3.2 mm high (0.5 x 0.125 in.) 

The compression test specimens were assembled to circular plates, with the 
bottom plate anchored to the snaker table and the top plate connected to the 
resonant mass. The assembly of the thirty 12.7 mm (0.5 in) high cylindrical 
specimens is shown in Figure 10. 

As with Che shear samples, the compression samples were each designed to provide 
similar stiffnesses. Appendix A discusses how the shape factors for the 
different height cylinders lead to selection of the above numerlca’ combinations. 

F. instrumentation 

The measurement requirements for the experimental Investigation of Che elastomer 
dynamic properties were as follows: 

1. Displacement measurement of elastomer support plate attached to the 
vibration table, relative to ground; 

2. Displacement measurement of elastomer support plate attached to Che 
resonant mass, relative to ground; 
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Fig. 10 Test Assembly of Thirty Elastomer Compression Specimens, 1.27 cm Diamet 
by 1.27 cm High (0.5 x 0.5 in.) 


3. Phase angle measurement between displacement measurements (1) and 
(2) above; 

4. Displacement measurement between two elastomer support plates, 
relative amplitude across the elastomer; 

5. Vibration frequency; 

6. Temperature of elastomer support plates and ambient air temperature 
in test specimen cavity; 

7. Vibrational power dissipation across the elastomer test sample. 

For convenience, acceleration was measured in (1) and (2) instead of displace- 
ment, but since the motions were sinusoidal the displacements were directly 
derivable from the accelerations. The very low level of motion of the shaker 
body (from which the shaker table is isolated by air springs with a very low 
resonant frequency) was also measured by accelerometers. 

The input acceleration (1) was measured by an accelerometer located on the 
support base for the elastomer test samples . Output acceleration on the 
resonant mass (2) was determined from three accelerometers mounted directly 
to the top plate attached to the elastomer specimens. The use of three 
accelerometers, spaced 90 degrees apart, with their output signals displayed 
on one oscilloscope screen, permitted immediate detection of nonaxial motions of 
the resonant mass. When nonaxial motions of the resonant mass occur, they 
manifest themselves either as amplitude or phase angle variations among the 
three signals, or as a combination of both (testing experience revealed that 
only the compression specimen having three cylindrical specimens was susceptible 
to nonaxial motions at certain test frequencies). 

The displacement measurement (4) between the two elastomer support plates 
(relative amplitude across the elastomer) was accomplished with a noncontacting 
capacitance-type sensor. 

Chromel-Alumel type thermocouples were used to measure the temperature at the 
midface location of two shear specimens in each test assembly. The thermo- 
couples were embedded in small holes in the aluminum holding plates at a 
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distance of approyimately 1.6 mm (1/lb in.) from the elastomer surfaces. 

The compression test samples contain thermocouples directly above and below 
the center point of the flat surfaces of three individual cylinders. Again, 
the thermocouples were separated by approximately 1.6 mm (1/16 in.) of metal 
from the elastomer surfaces. (See Figure lx for typical thermocouple instal 
lation in compression test sample.) 

Displacement and acceleration data, together with temperatures from individual 
thermocouples were displayed visually in digital form for monitoring purposes. 

As a permanent record, this data was also recorded on a digital printer, with one 
thermocouple reading automatically printed with each amplitude data point. The 
two types of data were cycled independently by their own scanner-switches, and thi s 
a given vibration sensor output did not always appear on the same line as a given 
thermocouple output. This is illustrated in Figure 12. 

Amplitude signals, either acceleration or displacement, were sequentially 
switched by an analog scanner into a two-channel tracking filter, which provided 
a visual readout of vibration frequency and of two filtered amplitude signals. 

The base amplitude signal (acceleration) was fed at all times into one of the 
two channels of the tracking filter and from there into a phase meter where it 
served as a reference signal for the measurement of the phase angle of the 
resonant mass amplitude signal (one of the accelerometer signals). The d-c 
values proportional to phase angle and amplitudes from the phase meter and 
the tracking filter were then converted into digital forr in digital voltmeters 
and printed on paper by a 21-column printer at a rate of approximately two lines 
per second (manually activated). Each line additionally contains the vibra- 
tional test frequency and the temperature indication from one thermocouple. 

(see Figure 12 for a typical data printout, Figure 13 for a photograph of 
the instrumentation console and Figure 14 for a schematic of the data acquisi- 
tion system) •• 

To monitor the power dissipation in the test sample during testing, a power 
meter was specially designed and built. This meter utilizes available 
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Fig* 14 Schematic of Data Acquisition System for Measurement 
of Elastomer Dynamic Properties 












acceleration and displacement signal inputs. The three (unfiltered) acceler- 
ometer signals from the top plate of the specimen holder (resonant mass signals) 
are averaged, to give a signal proportional to the net force on the resonant 
mass. This signal is therefore almost equal to that portion of the force 
acting through the elastomer sample which is due to dynamic displacement 
(differing only b} the small guide spoke and air cylinder forces due to their 
dynamic displacement.) The displacement signal which indicates relative motion 
across the sample is differentiated to give a signal proportional to rate of 
deformation (relative velocity). The force and velocity signals are multiplied 
together to give the Instantaneous power. The instantaneous power consists j 

of a d-c (average) component and a sinusoidally varying component and, by 
filtering out the sinusoidally varying component, the average power dissipa- 
tion is obtained. I 

W 

During test, the average power dissipation was held to a value which, by heat | 

transfer analysis, gave a temperature rise in the specimen of less than 

1°C. 

G. Testing Procedure 

One of the major features of the current test program was the performance of 
the scheduled tests at near-resonance conditions. Since phase angle between 
base excitation and resonant mass response is an accurate indicator of the amount 
of damping in the region of resonance, measurements should preferably be made in 
the phase angle range between approximately 15 and 163 degrees (see Appendix B) . 

This requires that the test frequency be approximately 0.9 to 1.5 times the 
critical frequency of the elastomer-resonant mass system. Coverage of the full 
test frequency range of 100 to 1000 Hz with sufficient test points therefore 
necessitates changes in the size of the resonant mass. Typically between seven 
and nine different mass arrangements ranging from 0.8 Kg (1.8 lb.) to 51 Kg (113 lb.) 
were found sufficient to cover the desired test range. The smallest mass consisted 
of the elastomer test sample top plate and the preload cylinder piston, with no 
additional weight attached. The next three combinations were obtained by attaching 
weights of up to 1.9 Kg (3.5 lb.) directly to the preload piston cylinder without 
the use of the spoke guide bearings. For the last five mass arrangements both 

guide bearings were used. 
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Once any one of the six elastomer test samples has been mounted on the shake 
table and the necessary instrumentation installed and connected, a typical 
test sequency proceeded as follows: 

1. A resonant mass was selected and installed. 

2. The elastomer sample was statically preloaded. (For zero preload 
cases, the static loading on the elastomer by the resonant mass 
was relieved through an appropriate pressurization of the top air 
cylinder. If the preload requirement exceeded the static deflection 
of the elastomer obtained from the resonant mass, pressurized air 
was admitted in the lower cylinder to produce the desired preload. 
Elastomer deflections were determined by readout from the capacitance 
probe installed in the lower base plate of the elastomer specimen 
holding fixture.) 

3. The elastomer test sample cavity was enclosed and the temperature 
control system given time to adjust ambient temperature to the 
desired value. 

4. With low vibration levels applied to the base of the elastomer 
holding fixture, frequency scans were conducted until the approximate 
resonant frequency of the system was found. Subsequently, an 

exact determination of the system resonant frequency was made by 
tracing, on an X-Y plotter, the d-c value proportional to the resonant 
mass acceleration amplitude as a function of the base excitation 
frequency. It may be noted here that, for a base-excited single 
degree of freedom spring-damper-mass system, resonance occurs at 
an angle smaller than 90 degrees. The deviation from 90 degrees is 
essentially determined by the amount of damping in the system as 
shown in Appendix B. 

5. While the predetermined vibrational power dissipation level in 
the elastomer test sample was maintained by adjustment of the 
shaker power input level, the vibration frequency was reduced 
until the phase angle decreases to a value between 15 and 25 
degrees. Provided none of the oscilloscope monitored accelera- 
tion and displacement signals showed signs of abnormalities (dis- 
tortions, or indications of nonaxial motion of the resonant mass), 2 g 



a data point was recorded. The recording of one data point 
consisted of either one or two complete printouts of all 
amplitudes (up to seven individual signals), associated 
phase angles and test frequency, and all temperatures from 
the thermocouples located in the test specimen holding plates. 

6. Stepwise increases in vibration frequency were imposed and 
between six and ten data points were recorded until the phase 
angle reached approximately 165 degrees. 

7. Shake table power was turned off and the preload on the 
elastomer sample was increased to the next higher value. 

8. Steps 4 through 7 were repeated until data at all desired pre- 
loads was obtained. 

9. Shaker table power was turned off, and all preload on the 
elastomer sample removed. The elastomer test sample was then given 
time (up to one hour) to recover its original, uncompressed 
height (as indicated by displacement probe readings). 

Tests, comprising steps 1 through 9, were then repeated with each of the 
remaining masses in turn, each mass giving a dynamic system with a different 
resonant frequency. 

The complete set of results obtained for each test specimen was reduced to 
yield stiffness and damping values over the range of frequency covered by 
the tests. The analysis upon which this data reduction was based is described 
in the following sections of this report. 

H. Test Rig Response 

With the forced vibration type of testing described in this section, the dynamic 
characteristics of the elastomer are basically determined by measuring the 
motions of the resonant mass and knowing the stiffness and damping characteristics 
of all elements of the test rig. It was, therefore, nece sary to analyze the 
dynamic response of the test rig with the objective of determining the character- 
istics of the various system components, viz., lower cylinder, upper cylinder and 
guide bearings. 
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A dynamic model of the cest rig is shown in Figure 15 where the system components 
are denoted by a spring or a combination of a spring and a dashpot depending 
on their physical characteristics. Both c> binders and the guide bearings have been 
designed such that their stiffnesses are negligibly small compared to the stiff- 
ness of the elastomer. Hence, very little precision is required in determining 
these stiffness values. However, damping introduced by the use of rubber seals 
in the cylinders needs careful attention since these damping coefficients may 
be quite comparable with that of the elastomer specimen. 

In order to determine the dynamic response of the test rig, each component was 
tested individually and its relevant characteristics were measured. When 
assembly does not introduce any appreciable changes, the complete dynamic re- 
sponse of the test rig is determinable by the characteristics of each of the 
system elements. The most likely elements, whose response may depend on the 
assembly, are perhaps the rubber seals in the cylinders since their damping 
coefficients may depend on their proper alignment to a certain degree. With 
the primary objective of precisely determining the dynamic response of the 
test rig, the elastomer specimen was replaced by a mechanical spring shown in 
Figure 16, whose stiffness and equivalent mass are measured experimentally. Along 
with the testing of each system element separately, the total response of the test 
rig with the mechanical spring was measured as a function of air pressures in the 
cylinders. Although the stiffness of the various elements may be determined 
by the experiments performed with each of the system elements, it was found 
that the damping of the cylinders was best determined by performing a least 
squares analysis of the measured response of the assembled test rig. 

The details of the analyses described above are presented in Appendix B. The 
results arc summarized below. 

Numerical Values for Dynamic Characteristics of Test Rig Elements 
Guide Bearing Stiffness » 3.675 x 10^ N/m (210 Ib/in) 

Lower Cylinder Damping c^ * 1.02 x 10^ N-sec/m (0.583 lb-sec/ in) 

Upper Cylinder Damping c^ * 2.45 N-sec/m (0.014 Ib-sec/in) 
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Lower Cylinder Stiffness » 1.642 x 10^ + 1.023 N/m (938 + 40.25 Pj,^|||r 

Upper Cylinder Stiffness = 9.0 x 10^ + 1.525 p^ N/m(515 + 60.00 Ib/in) 

2 0 

where p and p are pressures in N/m or Ib/in'^gage in the lower and upper cyl- 

inders respectively. Note that all the above values are independent of frequency. | 

As shown in Appendix B, when the elastomer specimen is in place^ and the phase m 

angle p and ratio of output to input amplitude, a, is measured, the elastomer I 

stiffness and damping (Ke,Ce) are calculated as follows; ^ 

2 2 I 

,, „ . (ku + 2ks)a(cos4)-a)+(jjCuaSin4>-maj (acos(p-a ) | 

Ke - -K. + 

X. 2 

a - 2acos <^ + 1 m 

( 1 ) I 

r . Cua(cos<{)-a)-Hna)otSln(t> (l- ^^ - 
ue - -c . + mw 

^ 1 

a - 2acoscJ) + 1 

All the data for damping of the elastomer is presented as the imaginary part ^ 

of the complex stiffness (a»Ce). | 

I 

I 

I 

I 

I 

I 

I 
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IV . EXPERIMENTAL RESULTS 


Dynamic response data for several elastomer samples was obtained over a range 
of excitation frequencies. The power dissipated in the elastomer was monitored 
within the present experimental procedures to establish an upper bound on the 
temperature distribution within the elastomer sample. The ambient temperature 
around the elastomer was maintained constant. Under such conditions, the data 
obtained determines the elastomer behavior at fixed temperature. 

A. Operating Conditions and Data Range 

Three compression and three shear samples as described earlier were tested. 

The preloads for the compression samples were set to 2.5 and 5 percent of the 
heiglit of the specimens. Under all operating conditions, the elastomer was 
never subjected to tension in the compression experiments. In the case of the 
shear samples, three different preload, viz., 0., 2.5, and 5 percent of the 
specimen thickness were used. The frequency range for cJ1 experiments was 100 
to 1000 Hz. The ambient temperature around the elastomer was maintained to 
32°C within l^C. Also, the total power dissipated in the lower air cylinder 
and the elastomer sample was held constant tc establish an upper bound on the 
temperature gradients in the elastomer sample. The conditions were adjusted 
such that the maximum temperature rise within the elastomer specimen was less 
than l^C as discussed below. 

B. Temperature Gradient in the Elastomer 

Since the elastomer sample and the lower air cylinder act in parallel, the 
total measured power dissipation is divided between the cylinder and the elasto- 
mer specimen, depending on their relative damping characteristics. This total 
power per unit of the elastomer volume was maintained constant at about 0.022 
watts/cm in all the experiements . For the purpose of computing an upper bound 
on the temperature rise from either end to the center of the specimen, it may 
be assumed that the entire power is dissipated in the elastomer and, hence, 
writing the energy equation gives 


where T is the temperature (^C), q is the dissipation (watts/cm^), is the 
thermal condictivity (watts/^C cm) and x is the coordinate as shown in Figure 17. 


Since the temperature at the two ends is held constant at say T » the above 

o 

equation may be integrated to determine the temperature T in the center of 

c 

the specimen. 



where a, height of the specimen is 1.27 cm (0.5 in.) for the worst case. Assuming 
the thermal conductivity to be 0.002 watts/^C cm [7] the temperature rise for 

3 

a dissipation of 0.022 watts/cm will be 


T - T 
c o 


0.022 

X 

8 X 0.002 


(0.50)^ 


= 0.344°C 


Thus, under all operating conditions for the entire experimental data, the 
temperature rise within the elastomer is less than l^C. Since the ambient 
temperature was monitored to be constant at 32^C to within l^C, it may be con- 
cluded that the measured elastomer response is under isothermal conditions. 

C. Reduction of Experimental Data to Complex Stiffness 

The experimental variables measured included the input and output acceleration, 
relative phase angle, relative amplitude of elastomer displacement and the 
dissipated power. When the damping coefficient for the lower cylinder is a con- 
stant, then the dissipation in the cylinder at any frequency is proportional 
to the relative displacement amplitude. Thus, when the total dissipation (cylinder 
plus elastomer sample) is held constant, the displacement amplitude is an 
indicator of elastoner dissipation. From the measurement of the input and out- 
put acceleration, and the observation that all motions are sinusoidal, the system 
transmissibility a, which is the ratio of output to input anplitudes, ma ^ be deter- 
mined. i^rom the measurement of the phase angle cp and the computed a, Zcuations 

(1) andTabJel may be used to compute the elastomer stiffness k and damping 

e 

coefficients c . 

e 
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Fig. 17 Schematic of Thermal Gradients 
in the Elastomer Specimen 


The physical characteristics of an elastomer are sometimes denoted in terms 
of a complex stiffness. Also since simulation of material behavior by mechani- 
cal viscoelastic models is generally performed in terms of complex stiffness 
or complex compliance, these quantities are defined below. 

Complex Stiffness K = K1 + iK2 

Complex Compliance H = HI - iH2 


Also, note the relationships 
H = 1/K 

K1 = k the elastomer stiffness coefficient 
e 

and 


K2 = wc^ where w is the frequency in radiand/sec. 


With the above relations, the complex stiffness is also written as 

K = k + i (u)c ) . 
e e 

From the computed k^ and c^, the real and imaginary parts Kl and K2 may be deter- 
mined and all the results are presented in terms of these parameters. 


D. Compression Tests 

The results of the compression experiments are shown in Figures 18 to 23. It 
is clearly seen that dependence of Kl or K2 on frequency is rather small under 
the thermal conditions ensured by the power monitoring system. Also, the de- 
pendence of the complex stiffness on preload is insignificant. The large 
scatter in Figures 22 and 23 is primarily due to the geometry of this particular 
sample. There are three pieces of elastomer between the plates in this case. 

It is known that in order to prevent any irregular motion of the resonant mass 
at least three support points are required. However, if only three supports are 
provided, then any misbehavior at any one of the supports will result in somewhat 
irregular motion. In fact, for this sample, substantial phase differences 
between the three output accelerometers were sometimes observed. 
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Fig. 18 Complex Stiffness of the Compression Sample With Thirty 
Specimens at 2-1/2 Percent Preload 
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Fig. 19 Complex Stiffness of the Compression Sample With Thirty 
Specimens at Five Percent Preload 


40 


■ 


mu 


5TIFFNEE5 <Lfc/IN> 

XD^ xa**" xo' 


EURSTtSHER SPECXMEN D.2i5XN«-Xa. 
MRTERXRU s PotybuksdiBn* 

LBRXXNB = CamprBsiian 
POWER = .ns »sakbs 

PREL.ORX sS.SQQM 
MERN TEMP = 32 , BO "C 


K = <KX> + i<K2> Ih/ia 


U1 




K1 



<HZ> 


K2 

Ztt 


Fig. 20 Complex Stiffness of the Cooqpression Sample With Ten 
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Fig. 22 Conq;>lex Stiffness of the Compression Sample With Three 
Specimens at 2'l/2 Percent Preload 
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Fig. 23 Complex Stiffness of the Compression Sample With Three 
Specimens at Five Percent Preload 
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Close examination of the results of Figures 18 and 23 show the values for K1 

to fellow a number of distinguishable lines whose slopes differ from the overall 

trends. These lines respond to. the different values of resonant mass, and it 

should be noted that each change In mass requires a reassembly of the test rig 

and basically a new set-up. These observations Indicating the values of K1 to 

be somewhat dependent on the resonant mass, primarily Imply that the inferred 

values of K1 ;o a certain extent are dependent of the natural frequency of the 

test rig. This dependence is less pronounced for K2. Although these observations 

are not completely understood at this time, it may be mentioned that, in general, 

the sensitivity of a resonant mass method to residual errors in the test procedure 

can be expected to vary with (— ) . 

0) 

n 

The variation of relative amplitude of deformation in the elastomer as a function 
of frequency (when the dissipation in the elastomer plus the lower cylinder is 
held constant) was found to be quite independent of preload. A typical varla- 
tioa of amplitude as a function of frequency is shown in Figure 24. The decrease 
in amplitude with Increase in frequency indicates that the power dissipated in 
the elastomer also decreases with Increases in frequency since K2 or c^u is 
fairly constant. Tt should, however, be noted that the computations of K1 and 
K2 are based on the assumption that the air cylinder damping coefficients are 
constant . 

E. Shear Tests 

Complex stiffness data for the shear experiments for the various specimens 
at different preloads are summarized in Figures 25 to 33. Tn general, the 
frequenc lependence is again small and also the effects of preloads appear 
to be insignificant. However, when compared with the compression 
data, it is seen that both the real and imaginary parts of the complex stiff- 
ness are somewhat higher in the case of shear tests. 

Again, due to the power constraints, the relative amplitude variation as a 
function of frequency is nearly identical for all the specimens. Typical 
results are shown in Figure 24. 
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Fig, 24 Variation of Relative Deformation Amplitude 

in the Compression and Shear Elastomer Sample 
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Fig. 25 Complex Stiffness of the Shear Sample With Tvtnty 
Specimens at no Preload 
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Fig. 27 Complex Stiffness of the Shear Sample With Twenty 
Specimens at Five Percent Preload 
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Fig. 28 Complex Stiffness of the Shear Sample With Eight 
Specimens at no Preload 
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Fig. 33 Complex Stiffness of the Shear Sample With Four 
Specimens at Five Percent Preload 
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F. Sunnnary and Correlation of Test Results 


Within the experimental range of frequencies the complex stiffness shows a 
fairly mild frequency dependence for all compression and shear experiments. 

Thus, for the purpose of estimating specimen stiffness, single values of the complex 
stiffness may be stated for the entire frequency range without any appreciable 
errors. The second and fifth columns in Table 1 give the mean values over the 
frequency i ige of the real and imaginary stiffness components for each sample 
under each preload condition. The similarity of these mean values at different pre- 
loads Indicates insensitivity to preload. The similarity in stiffness among 
elastomer samples which have different shape factors, but have been designed for 
similar stiffness (see design details in Appendix A), is also fairly similar. 

The weak frequency dependence of the data may be summarized in terms of the 
coefficients of a power law correlation function of the form 

Y = A w® 

where Y is one of the variables of interest (Kl, K2, or amplitude) or 



Thus the data falls on a straight line on log-log paper. The coefficients 
A and B for each data set are computed by conventional least squares regres- 
sion analysis and summarized in the third, fourth, sixth, and seventh columns 
of Table 1. These correlation terms may be used to determine the relaxation 
modulus of the material at constant temperature. Deformation amplitudes are 
also summarized in Table 1, both in terms of mean values over the frequency 
range in terms of the coefficie. of the same power law correlation 
functions. As a demonstration, the power loss correlation to the 1.27 cm 
(0.5 in.) high compression sample with 2.5 per cent preload is superimposed on 
the experimental data in Figure 34. 

A powerful technique for the correlation of elastomer data with marked frequency 
dependence is to fit mechanical viscoelastic models. For the present data, 
with its weak dependence on frequency, such techniques may not be employed veiy 
efficiently. However, as a demonstration of the approach, the behavior of the 





TABLE 1 

SUMMARY 01- ELASTOMER TEST DATA 




0/1 




MATERIAL • Polybut«dltn« 

TEMPERATURE • 32 dtg-C 

CUNPLEX STIR-NESS K • Kl ♦ 1K2 Xb/Xn 
CORRELATIUN FORM Y • A«OMEOA*i« 


SPECIMEN 





AMPLITUDE 


DETAILS 

MEAN 

A 

B 

MEAN 

A 8 

MEAN 

A 

8 

Ib/ln 



lb/ in 


In 



O.EOIN-30^C/2.5 

7.54IE404 

2.224E^04 

l.889E*0l 

6.347E403 

t.B73E404 -I.05IE-01 

I.05SE-03 

3.854E-0I 

-7.677E-OI 

0.501N-30/C/9.0 

7.a09E^04 

2.545E404 

1.458E-01 

S.S57E403 

I.507E404 -7.2I2E-02 

I.042E-03 

3.40IE-0) 

-7.529E-0I 

0.251N>IO/C/2.9 

S.042E^04 

2.S7IE404 

I.32SE-OI 

9.266E403 

2.809E404 >I,427E>0I 

3.806E-04 

2.040E-0I 

-$.082E*0t 

0.2S1N>IO/C/'S.O 

6.026E404 

3.577E404 

I.044E-0I 

9.447E403 

2.753E404 -l.saEB-OI 

3.924E-04 

I.7.4E-0I 

•7.889E>OI 

O.I2»lN-3/’C/'2.5 

7.662E404 

7.3I2E404 

6.368E-03 

I.I90E«04 

I.378E404 -2.004E-02 

2.901 E-04 

2.407E-02 

•6.2ISE-0I 

O.I2SIN*3/C/E.O 

0.692E«O4 

7.007E404 

3.227&-02 

l.295E^04 

1.854E404 *4.859E-02 

2.250E-04 

3.947E-02 

-6.854E-0I 

O.SItf>0 /S/0. 

I.70SE405 

S.330E^04 

1.S48E-0I 

1 . 720E^04 

3.487E404 -8.8S6E-02 

3.I86E>04 

2.478E-0I 

-8.473E-OI 

O.SIN-8 /S/2.S 

l.666E^0S 

S.003E^04 

1.S30E-0I 

l.737E^04 

3.I38E404 -7.S2IE-02 

3.498E-04 

4.989E-Ot 

-9.232E-0I 

O.SIN-a /S/S.O 

l.598E^05 

5.39IE404 

I.39IE-01 

1 .6S3E404 

2.678E404 -6.I42E-02 

3.462E>04 

2.I62E-OI 

-8.I9SE-01 

I.OlM-4 /S/0. 

1.453E405 

2.S69E«04 

2.I99E-0I 

l.d24EM>4 

2.S00E404 •6.9I2E-02 

3.497E-04 

2.844E-0I 

-8.903E-0I 

I.OIN-4 /S/2.S 

I.3SIE405 

2.SS0E404 

2.I3SE-OI 

l.63IE^04 

2.3S8E404 -4.8S2E-02 

3.799E-04 

2.902E-01 

-8.309E-0I 

1.0 lN-4 /S/5.0 

l.307E^0S 

2.2d3E404 

2.238E-OI 

1.563E404 

I.923E404 -2.647E-02 

3.820E-04 

2.639E-01 

-8.339E-OI 

0.2IN>20 /S/0. 

l.S2SE^0S 

3.706E404 

l.88IE*Ot 

l.689E^04 

5.700E404 -I.S48E-0I 

3.395E-04 

2.I68E-OI 

*a.220E-0l 

0.2tN«20 /S/2.5 

I.S49E«05 

3.33SE404 

1 .936E-0I 

1 .681 £404 

5.S67E404 -I.S76E-0I 

3.398E*04 

2.630E-0I 

-8.3S6E-OI 

0.2IN-20 /S/5.0 

t.452E40S 

3.689E404 

1.743E-OI 

1 .600E404 

4.833E4Q4 -I.406E-0I 

3.483E-04 

2.240E-OI 

-S.I96E*0I 


COMPLEX STIFFNESS K • Ml ♦ 1K2 Nt«/« 

CORRELATION FORM Y « A«OMEOA**B 


SPECIMEN 


************** K2 M*«*«9******* 

*#•**««**** 

AMPLITUDE 


DETAILS 

MEAN 

A 

8 

MEAN 

A 8 

MEAN 

A 

8 


Ndw/n 



N««/n 


m 



0.90IN-30/C/2.9 

I.321E407 

3. 8949406 

1.5899-01 

I.462E406 

3.280E406 -I.09IE-0I 

2.679E-D9 

9.790E-03 

-7.477E-0I 

0.50tH*30/C/5.0 

1 . 36BE407 

4.457E406 

I.498E-0I 

I.916E406 

2.640E406 -7.2I2E>02 

2.646E-09 

S.639E-03 

-7.529E-0J 

0.25IN-I0/C/2.9 

1 .408E407 

9.027E406 

1.329E-0I 

1 .623E406 

4.920E406 •I.427EH3I 

9.66BE-06 

9.IS0E-03 

•8.082E-0I 

0.25IN-I0/C/5.0 

I.406E407 

6.264E406 

1.044E-0I 

I.654E406 

4.822E404 -I.382E-OI 

9.96SE-06 

4.4S0E-03 

-7.889E-0J 

0.I25IN-3/C/2.9 

I.342E407 

I.28IE407 

6.370E-03 

2.084E406 

2.4I4E408 -2.004E-02 

6.394E-06 

4. II4E^ 

-6.2I6E-0I 

0.I25IN-3/C/9.0 

I.957E407 

I.227E407 

3.227E-02 

2.268E406 

J.247E406 -4.899E-02 

9.7I9E-06 

9.0I0E-04 

-6.894E-01 

0.91N-8 /S/0. 

3.I49E407 

9.339E406 

1.948E-OS 

3.0I2E406 

6.094E406 -d.SS6E-02 

8.092E-06 

6.294E-03 

-8.473E-OI 

0.9IN-8 /S/2.5 

2.9I8E407 

8.742E*0« 

1.930E-01 

3.042E406 

9.497E406 -7.92IE-02 

8.889E-06 

I.266E-02 

-9.232E-0I 

0.9lN-d /S/9.0 

2.799E407 

9.3699406 

I.39IE-0I 

2.899E406 

4.889E408 -8.I4IE-02 

8.794E-06 

9.491 E-03 

-S.t98E-0l 

I.OlM-4 /S/0. 

2.949E407 

4.499E406 

2.199E-01 

2.844E406 

4.903E406 -«.9I2E«02 

S.68IE-06 

7.223E-03 

-8.903E-0I 

I.OlN-4 /S/2.9 

2.366E407 

4.469E406 

2.I39E-01 

2.897E406 

4.I83E«04 -4.883EH92 

9.648E-06 

6.396E-03 

-S.309E-01 

1.0 IN-4 /S/9.0 

2.289E407 

3.963E406 

2.238E-01 

2.737E406 

3.3.9E408 -2.647EH)2 

9.703E-06 

6.692E-03 

-S.339E-0I 

0.2 IN-20 /S/0. 

2.946E407 

6.494E406 

I.0SIE-OI 

2.999E406 

9.984E408 -l.548E^I 

S.624E-06 

9.906E-03 

-S.220E-0I 

0.2IN-70 /S/2.9 

2.7I3E407 

9.B43E406 

1.936E-0I 

2.944E406 

1.028E407 -l.576E<«l 

8.631 E-06 

6.6S0E-03 

-S.3S6E-0I 

0.2IN-.0 /S/9.0 

2.942E407 

6.46IE406 

1.7439-01 

2.803E406 

0.4«4E«O« -I.40.E-0I 

9.399E-06 

9.691 E-03 

-S.I96E-0I 

N.8. Fr9qu«ncy 

UMEOA in tht 

corrtlftion functien 

if in rfd/Mc* 





A sample consists of several specimens of elastomer arranged in a prescribed 
configuration, 

"Specimen details" code given in the table is derived as fellows: 

(1) Compression Sample: (Specimen thickness) IN - No. of Specimens/C/ 

(X Preload). 

The diameter of all compression specimens is 1/2 in. 

(2) Shear Sample: (Specimen length) IN - No. of Specimens/S/ 

(I Preload) 

The specimen thickness and width normal to 
direction of deformation for all shear apecimens 
are respectively 1/8 in and 1.92 in. 
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data for the compression saiAple having 30 cylindrical elastomer specimens with 
1.27 cm (0.5 in.) diameter and 1.27 cm (0.5 in.) thick was correlated with a 
double Voigt model. The least squares technique described in Appendix D was 
used to determine the various coefficients of the model. Since the selected 
model is of the series type, the analysis is best performed in terms of complex 
complia,. is. The results are shown in Figure 35. It is seen that most of the 
stiffness is determined by the static spring, and that the stiffness and damping 
behavior of the model follows reasonably closely the mean level of the measured 
stiffness and damping values. Within this order of agreement between experi- 
mental data and model predictions, a large number of different types of models 
may be shown to fit the data. However, the use of mean values, or the simple 
correlation functions, is adequate in the present case. 
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Fig. 35 Typical Correlation of Experimental Data With a Double 
Voigt Model 
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V. DISCUSSION 


Llxamlnation of the experimental data indicates tliat, if the therma] state of the 

elastomer is lield constant, the complex stiffness of the material is relatively 

insensitive to the variation of frequency. It is expected that, if K2 or 

c^o) is a constant, then for constant power dissipation in the elastomer, the 

relative amplitude must also be frequency independent. However, in all the 

experiments, the power measurement s primarily used to establish an upper 

bound on the temperature gradient in the elastomei and, in fact, the measured 

power is a sum of elastomer dissipation and the dissipation in the lower 

cylinder. It may be shown [8] that, if the damping coefficient of the cylinder 

is a constant, and if the tvOtal work done per cycle is when the relative 

amplitude is x , then 
r 

2 2 
IT C.U)X^ + ITC U)X =* W 

i r e r 

The above equation may be rearranged as 

2 W/ ff 2x power 

\ “ K2 + * w(K2 + 

Since the imaginary part of the complex stiffness K2 * c^io was found to be fairly 
insensitive to the frequency, the amplitude will decrease as a function of 
frequency at constant power, as found in all the experiments. 

The uncertainties in the computed elastomer properties will be generally small 
when the nominal stiffness and damping values of the tesc samples are substan-- 
tially larger than those of the various elements of the test rig, as sfiown by 
Equations (B5) and Table Bl. A comparison of the elastomer coefficients shown 
in Table I with those of the test rig components clearly shows that the stiff- 
ness of any of the test rig elements is indeed small compared to the elastomer 
stiffness and, therefore, any errors will not significantly affect the computed 
stiffness o/ the elastomer. In case of the damping coefficient, however. 
Equations (B5) need further discussion. If the proper values of c^ and 
from Table Bl .re substituted in the second of the Equations (BS), the elastome.'* 
damping ef fici^ nt may be written as 
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" TTf;: 


o 




a (cos(|) - a) + a 


1 


c. + c 
I e 


u 

mo) 


in 03 


a - 2a cos4> + 1 



Since the experimental data was collected at conditions close to resonance 
where transmissibility a > > 1, and the residual stiffnesses are small compared 
to the elastomer stiffness, the above equation may be approximated as 


I e ^ sin4) __ __u 

mu) a nu‘) 


I 




1 


or 


mo) sin 


- 'u> 


Thus, an error in the estimation of c. and c will have an additive error in 
’ £ u 

or a linear frequency dependent error in K2. 

It is now clear that, if the nominal value of elastomer damping is substantially 
large in comparison with the sum of the cylinder damping coefficients, (c^ + c^), 
then any small uncertainties in the estimation of c^ + c^ will not appreciably 
eifect the elastomer damping coefficients c^* If the results summarized in 
Table 1 are reexamined, it is seen that K2 = c^u3 is fairly insensitive to 
frequency. The elastomer damping, therefore, decreases as a function of fre- 
quency and the greatest significance of errors in cylinder damping are to be 

expected at the highest frequencies. The lowest values of K2 for the compression 

6 3 

and shear samples from Table 1 are 1.462 x 10 N/m (8.347 x 10 Ib/in.) and 
2.737 x 10^ N/m (1*563 x 10^ Ib/in.) respectively. At a frequency of 1000 Hz 

2 

these values result in elastomer damping coefficients, c , of about 2.33 x 10 N-sec/m 

2 ^ 
and 4.36 x 10 N-sec/m (1.33 and 2.49 Ib-sec/in) respectively for the compression 

2 

and shear samples. The total residual damping (c + c ) is about 1.05 x 10 N-sec/m 

At U 

(0.6 Ib-sec/in.). Thus, uncertainties in the estimated residual damping may 
become significant in the extreme high frequency region. Under such conditions, 
methods both to reduce and more accurately dete' mine these residual damping 
coefficients should be considered in relation to future test rig modifications. 
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The effects of preloads, if any, on the dynamic response of the elastomer are 
not clearly noticeable within the present data. Of course the maximum nominal 
strain due to the preload was no more than 5 percent. At least for the com- 
pression tests it might be expected that an increase in preload results in 
increased stiffness. This is indeed observed to some degree as shown in 
Table 1, however, the present data does not warrant such a general conclusion. 
In fact, for the shear samples a reversed effect is consistently observed. 

Since vibration amplitude was not varied as an independent parameter, no con- 
clusions can be drawn as to the influence of amplitude. It is, ho^/ever, a 
premise of this work that, as long as the temperature of the elastomer remains 
constant, the influence of amplitude will be '^mall. In other words dynamic 
characteristics measured under conditions which both limit the temperature rise 
and control the environmental temperature should be the same whether the amp- 
litude is constant or the power dissipated is constant. Further investigation 
of the influence of amplitude will be necessary to consolidate this point. 

Relations between the complex stiffness and some basic physical properties 
such as relaxation modulus could be mathematically expressed [9] in terms of 
Fourier relations. However, a much more rigorous analysis of the experimental 
data will be required in order to obtain the correlation functions describing 
the complex stiffness in the entire frequency range of zero to infinity. This 
clearly warrants more experiments and data analysis. It is anticipated that 
future investigations will accomplish this objective. 

As described in Appendix A, the test samples were designed with the primary 
objective of Investigating the elastomer behavior in the frequency range of 
100 to 1000 Hz with varying shape factors. This objective has clearly been 
achieved in the present investigation. Furthermore, the results summarized 
in Table 1 show that variations in the dynamic stiffness within the three 
compression or the three shear samples are small. However, the dynamic 
stiffness of the shear samples is about twice that of the compression samples. 
If the static stiffness ratio of shear to compression samples is computed 


from the simple relations given in Appendix A, it is found that this ratio 
varies from 1.8 to 2.27. Thus, the observed difference in dynamic stiffness 
of the shear and compression samples is primarily due to the differences in 
the static stiffness for the selected geometrical configurations. The 
absolute values of the dynamic stiffness, as might be expected, are substantial! 
larger than those of the static stiffness (as defined in Appendix A). This 
is in agreement with most of the published investigations, e.g., Meyer and 
Sommer [5] or Cardillo [4]. Although no general conclusion with regard to 
the influence of shape factors on the dynamic properties of an elastomer 
sample may be presently justified, it is very likely that shape factors 
influence the static stiffness as described in Appendix A. In other words, 
samples with varying shape factors but similar static stiffness will have 
closely identical dynamic stiffness at any frequency, as demonstrated in 
this investigation. 


VI. CONCLUSIONS 


A BTR elastomer (polybutadiene) has been tested under shear and compressive 
loading, at 32°C, for constant power dissipation, over the frequency range 
100 to 1000 Hz. Displacement amplitudes covered 2.54 x 10*^ to 7.6 x 10”\ 

(strains of 0.0008 to 0.024). From the test results, for these conditions, 
the following conclusions may be drawn. 

1. The procedure for monitoring and controlling power dissipation and 
environmental temperature was successful in allowing measurement of 
sample properties at a temperature constant with + 1°C. 

2. The dynamic properties in terms of complex stiffness are only weakly 
influenced by frequency. 

3. A simple power law relationship may be used to represent the measured 
frequency dependence. 

4. The dynamic properties in terms of complex stiffness are not significantly 
influenced by preload. 

5. The use of published shape factors to lesign different shaped elastomer 
samples for equal stiffness resulted -*n samples whose dynamic stiffness 
and damping values were similar. 
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VII. RECOMMENDATIONS FOR FUTURE RESEARCH 


On the basis of the work reported and the overall objectives of the Elastomer 
Technology Program, the following recommendations for future investigations 
are made. 


1. Utilize the reported stiffness and damping properties for elastomer 
specimens of particular geometry to determine generally applicable 
material properties - 

2. Apply these material properties in the prediction of stiffness and 
damping properties of the same material, but with different geometries. 

3. In particular, predict the stiffness and damping of a cylindrical 
cartridge (such as might be used to flexibly support a rolling element 
bearing) . 

4. Test a cylindrical cartridge to determine its dynamic stiffness and 
damping properties and compare these measurements with corresponding 
predictions. 

5. Design a cylindrical cartridge damper for controlling vibration of 
a specific rotor-bearing system and test the performance of the 
cartridge in this application. 

6. By further controlled-temperature tests on elastomer specimens, deter- 
mine what influence temperature has on the material properties of 

the BTR elastomer tested in the presently reported iivestigation. 

7. Develop a model for prediction of elastomer behavior which accounts 
for temperature dependence of material properties. 

8. As identified in this report, uncertainties in the damping of the 
lower cylinder could be significant at the higher test frequencies. 
Reduction of this uncertainty should be factored into future rig 
modifications by such means as use of a metallic foil to ^.iplace 
the present rubber seals. 


APPENDIX A 


TEST SAMPLE DESIGN 

The frequency range over which the elastomer samples are required to be 
tested was a dominating factor in the design of the samples since accept* 
able data cannot be collected at frequencies far from the resonance. This 
may be seen by examination of a typical response plot for a base-excitation 
system, as described in Appendix B. For frequencies greater than three times the 
natural frequency, the system transmissibility approaches 1.0 and, depending 
on the amount of damping in the system, the phase angle between the input and 
output displacements tends to 180 degrees. Based on these facts, it was de- 
cided that no data should be collected at frequencies higher than three times the 
natural frequency. Since from a practical standpoint the upper bound on the 
vibration frequency is on the order of 1000 Hz, it was decided that the 
test set-up should provide a resonance frequency of at least 350 Hz with the 
smallest mass, which was approximately 0.9 Kg (2 lb.) in the present investi- 
gation. These constraints provided a stiffness which the test sample must 
have in order to meet the experimental requirements. Both the compression 
and shear samples were designed based on this criterion. 

1. Compression Sample 

The compression sample is an ensemble of several cylindrical elastomer speci- 
mens, as described in Section III. The size and number of these specimens was 
determined by the above design criterion. If d is the diameter and "a’* is 
the height of any cylindrical elastomer specimen (as shown in Figure Al), 
and if a compressive load is applied along the axis of the specimen, then the 
nominal compressive stress is given by Payne and Scott [10] 


G 


- X 


I X^ 


(Al) 


where G = Shear modulus of the material 

X == Ratio of strained to unstrained length 
S = Shape factor 
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The strain, c, is simply (1 - X) and therefore Equation (Al) may also be 
written as 


o = G 


[a - e)-2 


1 + e 


(A2) 


Since for most of the experiments in the present investigation the small 
strain assumption will be valid, Equation (A2) may be expanded in powers of 
e, Hence 


0 = GS[3g + + ••.] 


(A3) 


If higher order terms are neglected in the above equation, then 


0 = 3GSe 


(A4) 


It should be noted that Equation (A4) is basically Hooke's Law, since the 
Young's modulus for most rubbers, E = 3G. 

The static stiffness of the elastomer sample is defined by load per unit 
deflection, i.e.. 


oAn 

St ea 


(A5) 


where A is the cross sectional area of each elastomer specimen and n is the 
number of specimens in the ensemble. 

Combining (A4) and (A5) results in 


k = 
s 


3GSAn 

a 


(A6) 


Equation (A6) gives the static stiffness. However, in dynamic testing the 
natural frequency of the system will be primarily determined by the dynamic 
stiffness. Meyer and Sommer [5] have produced test results comparing the 
static and dynamic stiffnesses of several Polybutadiene compositions. For 
a composition having a Shore A hardness of 71, it is shown that the dynamic 
stiffness is three times larger than the static stiffness. Since no 
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compositional details on the material used in this investigation are available 
except that it has a Shore A hardness of 70, the above findings is used to 
compute the dynamic stiffness, k. 


k = 3k 


9GSAn 


The natural frequency for a mass m is given by 


2 k 9GSAn 


where the cross sectional area is represented in terms of the diameter d. 


The minimum mass for the present system is 0.9 Kg (2 lb) and the desired 
is 350 Hz. Also, as a practical value, 1.27 cm (0.5 in.) may be arbitrarily 
selected for the diameter. Substitution of these values in Equation (A8) 
results in 

^ 2(350 X 2tt)^ 

a 9fr G X 0.25 x 386 


1.415 X 10 
G 


Note that the shear modulus G is in Ib/in and the height of the specimen "a" 
is in inches in Equation (A9). 


It is now only necessary to estimate the shear modulus and determine the 
shape factor S. Gent [11] has shown an approximate inter-relationship 
between Shore A hardness, s, and the Young's modulus of elasticity, E. 
Assuming G*E/3 this relationship is given by 


I. 


G 


(AlO) 


14.22 (56 + 7.66s) 

3 X 2.67 X 0.0515 x (254 - 2.54s) 


X 6.894 X 10 N/m 


Equation (AlO) is plotted in Figure A2, taken from Gent. Using thie rela- 
tionship, the shear modulus for a 70 durometer or Shore A hardness is estimated 
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as 1.86 X 10 N/m (270 Ib/in ). Substitution of this in Equation (A9) gives 

= 52.5 (All) 


The shape factor, S, is defined by Payne and Scott [10] by a function 


S 


1 + B 



2 


6 2 

where B is a constant which depends on the shear modulus. For G = 1.86 x 10 N/m 
2 

(270 Ib/ln ), B is estimated as 0.063. Thus, 


S 


1 + 0.063 



a i 


With the objective of studying the effects of shape factor on the dynamic 
properties, three different thicknesses of 12.7, 6.35, and 3.675 mm (0.5, 

0.25, and 0.125 in.) were selected to give shape factors of approximately 1, 

1.25 and 2. For these thicknesses and the shape factors, the number of specimens 
required in each ensemble were determined to be respectively 25, 10, and 3 for 
the three samples. For the first sample, the number n was in fact changed to 
30 primarily because of pratlcal reasons for geometrical symmetry of the sample. 
In any case this will increase the natural frequency and is, »-herefore, perfectly 
acceptable. 


2 . Shear Sample 

The design of the shear samples is somewhat simpler than that of the compression 
samples since no shape factors are associated with shear deformation. However, 
as shown in Figure A3, when a piece of rubber is shear-strained in a conven- 
tional manner both shear and bending occur. It has been shown [10] that the 
applied shear force, F, results in a deflection, 6, given by the expression 
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(A13) 


- Ik 

GA 


where b is the specimen thickness as shown in Figure A3, G is the shear modulus, 
p is the radius of gyration of the cross section about the neutral '^xis of 
bending and A is the cross sectional area. The bending effect is lndi..ated 
by the second term in the brackets. 

2 2 

For the rectangular sample shown in Figure A3, A “ wz and p z /12. By 
substituting these quantities in Equation (A13), the sample static stiffness 
in shear is expressed as 


I- 36d^ J 


(A14) 


2 2 

Bending will be ne^^liglble when b /36p < < 1, and for this simple case 


(A15) 


As in the case of the compression specimens, assuming that the dynamic stiff- 
ness, k, is three times the static stiffness, the sesonant frequency o) for 

n 

a mass m and a sample consisting of n identical specinens in parallel is 
given by 


3Gwzn 

mb 


(A16) 


Now if a resonant frequency of 350 Hz is desired for the minimum possible 
mass of 0.9 Kg (2 lbs). Equation (A16) gives 


2 X (350 X 2ny 
' 3 X 386 X 270 


- 31 


(A17) 


If the ,\ctual value of wzn/b is preater than 31. then the resonant frequency 
will be somewhat greater than 350 Hz (as given by Equation A16) . Such a 
design will be perfectly acceptable as long as .e lowest desired resonant 
frequency is attainable with the maximum avail •^hle mass. With these design guide 
lines, it was decided that a satisfactory shear sample may consist of four 
elastomer specimens with the following dimensions 

Sample No, 1: n = 4; w 5.08 cm(2 in.); z ■ 2.54 cm(1.0 in.) 

It may be noted that the above configuration gives wzn/b * 64 and therefore 

the resonant frequency corresponding f * the mass cf 2 lb will be 35C 64/31 

or 502 Hz. If the lowest desired frequency is 100 Hz, then the re/ 'red mass 

will be approximately 51 lb, which is well within the experimental limitations 

on maximum mass. Also, bending effects have been neglected in the above 

2 2 

design, which is quite justifiable since the bending term (b /36p“) in this 
case is 0.0052, which is indeed quite small compared to unity. 

In order to allow for some bending effects vlthin the experimental liinl tat ions, 
two additional configurations were selected such that (wzn/b) is held constant 
for all three samples to meet the stiffness requirements. The selected con- 


figurations are 









Sample No. 2: 

n » 

8; w * 5.08 cm (2.0 

in. ) ; 

b » 0.318 cm 

(0, 

,125 

in.) 


2 ® 

1.27 cm (0.5 

in.) 






Sample No. 3; 

n * 

20; w - 5.08 

cm(2.0 

in.); 

t " 0,318 cm 

(0, 

,125 

in.) 


2 * 

0.'08 '.m (0.2 

in. ) 







2 2 

The magnitude of the bending term (b /36p ) is estimated to be 0.0208 and 0.13 
for Samples No. 2 and No. 3, respectively. Thus the shear specimens have seme 
variation in the possible bending effects. 
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APPENDIX B 


DYNAMIC CHARACIERISTICS OF TEST RIG CC»dPONENTS 


Guide Bearings 

The two guide bearings, consisting of spoke assemblies, are identical with re- 
gard to their geometry and dynamic response^ Therefore, the experiment 
discussed below applies to both spoke assemblies- 

Since the only mode of dissipation in the spoke assemblies is internal friction 
in the metal and at the joints, the associated damping will be negjigibly 
small. It was found thac within the desired precision on the test rig response, 
the damping of the spokes may be neglected. In order to measure the stiffness 
a simple experiment was performed. The outside ring of the spoke assembly was 
clamped on the shaker table, as shown schematically in Figure Bl. The natural 
frequency of the system was measured, and k. . wing the mass of the hub* (see Sec- 
tion III), the stiffness of the spokes was computed in a straight forward manner. 
The duced dynamic model and the system c>iaracteristics are also outlined in 

Figu’-o Bl. It was found that the stiffness of each spoke assembly, k , is 

4 ® 

3.67 X 10 N/m (210 Ib/in). 

B. Mechanical Spring 

Any damping in the mechanical spring is once again negligible. However, due 

to th, size of the spring (see Figure 16), its equivalent mass must be 

determined experimentally. The spring was clamped on the shaker table and 

the natural frequency w of the system was measures. The dynamic model, in 

m 

this case, will be similar to the one sh^wn in Figure Bl. In the second experi- 
ment, a small mass, A was attached to the top of the spring and the natural fre- 

m 

quency was again measured. From these measurements, the computation of equi- 


valent spring mass, 

m 

and the stiffness, k is straightforward 

m 

“m J 

nr 

f 

and 

/k 

u’ - / *“ - 

m 

ro \/ m + Am 

V 

m 


y m 
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*The equivclent mass of the spokes is negligible compared to the mass of the hub* 




Dynamic Model of the Guide Bearing 



Thus 


and 


2 

o) m 

m m 


(Bl) 


Am 




Am 0) 

m 



It was found that the equivalent mass m = 0.877 kg (1.93 lbs) and the 

m 

corresponding stiffness = 2.96 x 10^ N/m (1695 Ib/in). 

C. Lower Cylinder with Free Pist m 

It is possible to mount the lower cylinder with the mechanical spring 
directly on the shaker table and obtain a base excitation system, shown 
schematically in Figure B2. The analysis for such a system is quite 
simple [12], and if the output motion is measured as a function of the 
input table excitation, the stiffness and damping of the cylinder can be 
calculated. It is assumed that the stiffness of the cylinder is a function 
of the air pressure only. 


The analysis for the base excitation system of Figure B2 is outlined in 
Appendix C, and it is shown that if the system transmissibility , a (defined 
by the ratio of output to input amplitude) and the phase angle between 
the output and input are measured, the values stiffness k^ and c^ may be 
computed. Within the scope of the present investigation, it is conserva- 
tively assumed that k^ and are frequency independent, as discussed in the 
Discussion section. Under this assumption, it is shown in Appendix C that 
only the measurements of a and the frequency at resonance are necessary in 
order to determine k^ and c^. In fact, the phase angle (\> could be computed 
by knowing k^ and from the measurement of a at resonance, and the reso- 
nant frequency. Thus, the resonant frequency and values of a were measured 
as a function of the air pressure in the cylinder. The phase angle 4> at 
resonance was also measured for the purpose of comparing it with the 
computed value. The amplitudes were low enough such that fluctuations of 
pressure in the cylinder were less than 1% under the adiabatic assumption 
for the highest amplitudes imposed. 
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Dynamic Model of the Lower Cylinder Mounted on 
the Shr>ker Table '"ith the Mechanical Spring 


With the above measurements and the analysis described in Appendix C, 
and were computed, and finally, using the relations shown in Figure B2, k 
and c ^ were determined. The results are shown as a function of air pressure 
in the cylinder in Figure B3. The computed and measured phase angles at 
resonance are also compared in this plot, where the solid curve represents 
the computed results and the experimental data is shown by the discrete 
points. It is seen that within the precision limits of phase angle measure- 
ment the agreement is fairly good. The cylinder stiffness is very well 
approximated by a linear function of pressure and the variation in the 
damping is rather small. However, it should be noted that the motion of 
the piston is free in this experiment while it is constrained in the assembled 
test rig. Such a constraint may influence the cylinder characteristics, 
since it certainly influences alignment of the seals. 

In order to put some constraint on the piston motion, the lower guide bearing 
was clamped to the piston and the experiment was repeated. The analysis of 
this setup, and of a similar setup for the upper cylinder, are special cases 
of the general response of the test rig and hence this response is formulated 
before discussing the above two special cases. 

b. Generalized Test Rig Response 

The dynamic model of the test rig shown in Figure 15 is reducible to the 
form shown in Figure B4 by combining the various elements appearing in 
parallel. Now the equation of motion of the resonant mass will be incependent 
of the shaker air spring and it may be written as 

mx2 + c^(x2“X^) + c^x^ + k^Cx^-x^) + k2X^ = 0 (B2) 

It is found experimentally that all motions are sinusoidal and hence solution 
of the above equation of motion becomes straightforward 

x. ® sin ojt 

(B3) 

^2 ~ ^2 ” ♦) 
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Combining Equa^*ions (B2) and (B3) and separating the coefficients of sin wt 
and cos wt terms gives 


^2 

- sincp + — stncp + 

raco 


1 2 

coscp ^ + — 2 “ coscp = 

mu) 


and 


(154) 


"1 I 1 2 ^ •^2 

— cosf - a M- ‘^os'p — 2 — 2 

' f tno3 mo) 


= 0 


where 


or = A^/A^. 


Equations ^B4) form the two algebraic equations which may oe solved for any two 
unknowns. For the purpose of reducing the elastomer data, the unknowns will 
bt' and while in the experiments with the upper cylinder and the mechani- 
cal spring assembly, the unknowns are and c^. Also for calculating the test 
rig response the solutions are required for or and cp* Equations (B4) are, there- 
lore, rearranged in three different forms. 


mcD 


'2 2 2 
^ j a (cossp - u) + ~ or sincp - or cos^p + or 


n\o) 


or - 2 or cos^p + 1 


i 

mo) 


^2 ^2 
— or (cosep - or) + cr sincp 1 - — r 
mu) ^ ^2 

mco 

2 

cr - 2 cr cos^ + I 


mco 


1 + (££i2 . 1) . ll 
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(B7) 


The expressions of and in terms of 

shown in Figure 15 are summarized in Table B1 for 


the various coefficients 
the different test setups. 


It is clear that the experiment with the lower cylinder assembled with the 
guide bearing and the mechanical spring, is a special case of the above 
formulation. Likewise, the experiment with upper piston connected directly 
to the mechanical springs and the guide bearings forms another special case. 
Since, in the case of both of these experiments, the stiffness and damping 
coefficients are assumed to be frequency independent, the formulation could 
be expressed in terms of the resonance conditions. This will be very similar 
to the analysis presented for the base excitation system in Appendix C. 


The system transmissibility ot is written in terras of the nondimensional 
quantities. 


n - ^ . 

W = — : u, 

JO * 1 
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mao 


rau) 


moo 


raoo 


0) 


or 


2 ^ o2 2 


^ * C? (Vj - 2(uj + Uj) + (u^ + Uj)' 
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TABLE B1 


REDUCED STIFFNESS AND DAMPING COEFFICIENTS 


Test Setup 

'll 



h 

Mechanical Spring with 
All System Components 

k, + k 
Z m 

^Z 

k 

u 

+ 2k 

s 

Elastomer Specimen with 
All System Components 

k. + k 
Z e 


k 

u 

+ 2k 

s 

Mechanical Spring with 
Lower Cylinder and 
Guide Bearing 

k + k, 
in Z 

^z 


k 

s 

Mechanical Spring with 
All System Components 
Except the Lower Cylinder 

k 

m 

0. 

k 

u 

+ 2k 

s 


Since at resonance a is a maximum, the second required equation is obtained 
by differentiating Equation (B8) with respect to U. 



or 


vj + 2 




(Vi + v^) - 2(u^ + 


U2> - 






(B9) 


Equations (B8) and (B9) are solved for either u^ and or u^ and Once 

these solutions are obtained, the phase angle at resonance is obtained by the 
second of Equations (B7), It is clear that the set of algebraic equations 
are nonlinear and, therefore, the conventional Newton-Raphson iteration 
method is used to obtain a final solution* The two equations are expressed 
in the following functional forms 


f 


f 


I 


2 


+ 2 uj ^ ■ '^1 


(BIO) 


Note that, when Equations (B8) and (B9) are satisfied, f * f^ = 0. 

The recurrence formulae for the computation of v^ or are given as 


u. - ^4 

^k+l k ‘k 




where 
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Once the above solutions are obtained, the phase angle at resonance is ex- 
pressed as 


tan;p 


(U2 '^2^ - n 

if (U, . vj - Vjvp - Uj (Uj + up 


(B12) 


rlius, if the stiffness and damping coefficients are independent of frequency, 
ii is only necessary to measure the system transmissibility and the frequency 
M 1 - soiKiiicc in order to compute the unknown coefficients. The special cases 
ot liu lower aiul upper cylinder with constraine d pis ton motion may now be con- 

s i<i . ' t 0 d . 


E. Lower Cylinder with Constrained Piston 

Some constraint on the motion cf the piston in tfa lower cylinder is provided 
by connecting it to the lower guide bearing. The dynamic model of this 
system is shown schematically in Figure B5. For a prescribed amplitude of 
table excitation, the frequency was varied until resonance was observed. 

The output motion of the resonant mass was then noted, and thus the trans- 
missibility a and the frequency at resonance were determined as functions 
of the air pressure in the lower cylinder. The amplitudes were low enough 
that fluctuations in the air pressure were than 1% assuming adiabatic 

expansion for the largest amplitudes impor ictions (Bll) with i - 1, 

were solved for Uj and v, and thus cyl < tiffnrss and damping 

coefficients were determined. 

The results are plotted as functions of preasu xgt re B6. Although 

the stiffness again varies linearly as a function of pressure, the absolute 
values are somewhat higher than those obtained with the free piston motion 
in Figure B3. A possible reason for this effect is the change in cylinder 
volume between the two tests [12J. The linear relationship shown by the 
dashed curve was used in calculating the lower cylinder stiffness for the 
purpose of reducing the elastomer data. 
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Variations in the damping coefficients are rather irregular, and assuming a 
constant value as shown by the dashed curve seems very reasonable. It 
should be noted that constraining the piston motion does Increase the damping 
to an appreciable extent. For reducing the elastomer data, cylinder damping 
was therefore determined by performing a least-squares analysis of the 
response of the completely assembled test rig, except that the mechanical 
spring was used instead of che elastomer specimen. Details of this approach 
are discussed later in this Appendix. 

The differences between the computed and measured phase angles are cf the 
same general order as shown earlier in Figure B3. 

F. Upper Cylinder with Constrained Piston 

Similar to the lower cylinder, an experiment was set-up with the upper 
cylinder mounted in place. The piston was connected through both guide 
bearings directly to the meachanical spring mounted on the shaker table. 

Thus, the lower cylinder was eliminated in this test. A dynamic model of 
the system is shown schematically in Figure B7. 

Since the range of pressure variation in the upper cylinder during the 
actual tests was much smaller than that for lower cylinder, the air pressure 

A 

was varied only up to 2.76 x 10^ N/m (4 pslg) in this experiment. The 
transmissibilitv a and the frequency at resonance were measured as functions 
of air pressure in the cylinder. Equations (Bll) with i«2, were solved in 
this case for U2 and v^ and thus the cylinder stiffness and damping were 
determined as functions of the pressure. 

The stiffness is again found to vary linearly with pressure as shown in 
Figure B8. The relationship shown by the. dashed curve was used during 
reduction of data from the elastomer tests. It may be noted that the 
stiffness of the upper cylinder is somewhat less than that of the lower 
cylinder. An explanation of this difference lies in the fact that in the 
assembled configuration, the volume in the upper cylinder is greater than 
that in the lower one. 
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The computed and measured phase angles In this experiment are In reasonable 
agreement and the order of deviations are the same as those shown earlier 
in Figure B3. 

For any finite air pressure, the damping vas again found to be fairly 
constant. However, the cylinder seems to provide significantly more damping 
at atmospheric pressure. After completing all the tests, the cylinders 
were disassembled to examine any differences between the mountings of the 
seals on the pistous. It was found that unlike the lower cylinder, the 
seal was not cemented to the entire area of the piston in the upper cylinder. 

This geometrical difference will affect the alignment of the seals and, 
hence, may explain some differences in the damping coefficients. In any 
case, the damping coefficients required for reducing the elastomer data were 
determined from the completely assembled test rig response as discussed 
below. 

G. Test Rig Response with a Mechanical Spring 

All components of the test rig were properly assembled in this experimental 
set-up except that the elastomer specimen was, replaced by the mechanical 
spring. The only damping in the system so assembled is due to the seals in 
the upper and lower cylinders. Also, the geometrical configuration of the 
test rig was very close to the one used when obtaining the elastomer response 
data. Thus, the alignment of the seals was very close to actual conditions 
and the computation of the damping coefficients from the response measurement 
of this system becomes quite meaningful. 

The dynamic model for the complete test rig was shown earlier in Figures 
15 and B4. The reduced form shown in Figure BA is directly applicable to 
the present configuration with the relations of the various system elements 
defined in Table Bl. The value of the resonant mass, m, was determined by 
summing up all relevant masses and was found to be 6 Kg (13.39 lbs). Stiff- 
nesses of the spoke assemblies, the mechanical spring and both the cylinders 
have already been established in the preceding experiments. Thus, from 
the measurement of the test rig response the damping coefficients for the 
upper and lower cylinders may be calculated. 
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The experimental pro .edure consisted of measuring the transmissibility a, 
the frequency, and the phase angle at resonance as functions of pressure in 
the cylinder. Pressure in the lower cylinder was varied from 0 to 8.27 x 
10^ N/m^ (0 to 12 psig) at three different values 0, 1.33 x 10^ and 2.76 x 
10^ N/m'^ (0, 2, and A psig) of pressure in the upper cylinder and thus, 
large data matrices were obtained to cover the entire range of operation of 
the test rig. Of the three measurements, a has probably the best accuracy. 
Frequency measurements are generally good to within 1 Hz but since most of 
the data obtained in this experiment lies in a very narrow band close to 
the natural frequency of the mechanical spring, small frequency errors can 
be significant. Keeping all these points in view, the damping coefficients 
were determined by a least-squares analysis of the transmissibility data. 

It is clear that the problem at hand is non-linear and some type of iterative 
or marching technique is required. The analysis used is briefly outlined 
below. 

A data matrix for c and c. varying from 0 to 1 Ib-sec/in was formed and for 

u ^ 

each combination, a at resonance is computed by Equations (B8) and (B9). 

Symbolically, let p and p. denote the pressures in the upper and lower 

cylinders where 1 ^ i ^ n and 1 ^ j ^ n^^ as determined by the selected 

experimental data matrix. Also, let c. and c have the values c. and c„ 

Jt u *’k ^ 

respectively where 1 ^ k ^ n and also 1 ^ h ^ n. If the measured trans- 
missibility for pressures p^^ and po is , and the computed value under 
these conditions and with = Cjj^ and c^ ® Cy^ be then the total 

squared deviation is obtained by proper summation: 

n n 

u I 

^h “ Z I ^ijkh' 

i-1 j-1 

The entire matrix for and Is scanned for the minimum squared deviation 

and the corresponding values of c^ and c^ are taken to be the required values. 

It is true that the accuracies of c^ and c^ are determined by the Internal size 

used to divide the region of 0 to 1 in finite number of data points. Initially 

the interval size is taken as 0.1 and the required c and c are determined and 

^ X 
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then a new matrix with interval size 0,01 is set up around tiiese computed 
values and the previously obtained solution is refined by adding an 
additional digit. The procedure could be repeated for obtaining any more 
significant digits if necessary. 

Using the above scheme, the damping coefficients for the lower and upper 
cylinders were found to be respectively 102 and 2,45 N-sec/m (0,583 and 
0.014 Ib-sec/ln), These values were used in determining the elastomer 
characteristics from the experimental data. 




APPENDIX C 


RESPONSE OF A BASE-EXCITATION SYSTEM 


Since the stiffness and damping in the air cylinders have been assumed to be 
independent of frequency, at least within the scope of the present investiga- 
tion, it is necessary to analyze the test rig response to determine these 
coefficients of damping and stiffness. In the absence of the upper cylinder 
the test rig is represented by a conventional base excitation system. The 
response of such a system is very well known [12] and this Appendix briefly 
outlines the analysis applicable to such a system. 


The response of a simple base excitation system, as shown in Figure Cl, is 
given in terms of transmissibility a, and the phase angle cp, defined by the 
equations given in Reference [9]. 


a 


and 


tan cp 





(Cl) 


where: 


Ui 


CO 


n 

k 


c 

c 


excitation frequency (rad /sec) 

the natural frequency of the spring mass system (rad/sec) 
ni 

spring stiffness (N/m) 


resonant mass (N-sec /in.) 


damping ratio 


system damping coefficient (N-sec/m) 
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Fig. Cl 


Pchetiifitic of a Base Excitation System 



damping coefficient (N-aec/m) 
■ input amplitude (m) 

== output amplitude (m) 
cp * phase angle (degree) 


Thus, if the Input and output amplitudes and the phase angle, cp, are measured, 
the system stiffness and damping may be computed from the above equations. For 
simplicity, Eq. (Cl) may be rearranged for this computation. 


k 

2 

OKU 


and 


c 

mu> 


Of (g - cos (d) 

1 + a (a - 2 cos cp) 


- o; sin CO 

1 + g (g - 2 cos cp) 


(C2) 


Since the measurement of the phase angles may sometimes Involve large errors, 

the values of k, and c, as computed by Equation (C2) may result in some error. 

However, the measurements of A 2 and A^^ or g and the frequency u) are relatively 

precise. Hence, If k and c are derived from the amplitude ratio g and the 

tnax 

frequency Q at resonance, then the errors in the computed k and c will be small 
compared to those determined by Equation (C2) . Now if k and c are assumed to 
be Independent of frequency, then the two equations required tc> compute k and c 
become: 


max 



and 



0, at CD • fl. 
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The preceding two equations are easily solved for (*^) and a as a function of 


CD 


max 


Q. The phase angle at resonance Is then computed from the second of Equation (Cl) 
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sr - 1 
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These solutions are readily plotted In Figure C2 with an additional notation: 
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(C5) 


The computations of k and c with the measured cr . and flare now straightforward 

IDsX 

from Figure C2 . For example, for a ■ 2.5, the following values are noted: 

Q » 0.22 ; 0 * 0 . )*2 ; 0 ■> 34 ** 

and from the definitions of ^ , 0, and 6, it Is found that cp " 56°, ■ 0.958. 

4° 

For a mass equivalent to say 0.9 kg (2 lb), this gives k ■ 7 .66 x 10 N/m (437 
Ib/ln.) and c ■ 37.1 N-sec/m (0.212 Ib-sec/ln.). 

Once k, and c are determined, the system response Is directly computed, using 
Equation (Cl). This response may be compared with the one obtained experi- 
mentally to ascertain some confidence on the computed values of k and c. 







APPENDIX D 


MECHANICAL MODEL ANALOGIES OP ELASTOMER DYNAMIC RESPONSE 

During the course of this study, techniques for fitting equivalent mechanical 
models to elastomer stiffness and damping data were developed. This Appendix 
describes these techniques. Figure 35 of the ijport compares the prediction of 
one such model with the data to which the model was fitted. 

Mechanical models simulating viscoelastic behavior of elastomers may be obtained 
by combining several springs and dashpots in a prescribed fashion. The simplest 
arrangements are when only one spring and one dashpot are connected either in 
series or in parallel as shown in Figure Dl. The series arrangement is generally 
called a Maxwell element, while the parallel configuration is termed as a Voigt 
element. Sometimes the parallel arrangement is also referred to as a Kelvin 
element, since it was Introduced by Kelvin before Voigt. In general, for a solid 
material, a Maxwell model may be defined by a combination .veral Maxwell 
elements connected in parallel with a spring and a dash .ikewise, a general 

Voigt model may consist of several Voigt elements conne ' i.n series with a 
spring. These models are shown schematically in Figure Although several 

other combinations are possible, Lazan [l3l has shown that a large number of 
arrangements are reducible to either a Maxwell or a Voigt type of model. The 
analysis for both these models is quite simple. In fact, if the Maxwell model 
is analyzed in terms of complex stiffness, the analysis will apply directly to 
the Voigt model when the stiffness function is replaced by compliance. 

1. Maxwell Model 

Since the force in a spring and the dashpot in a Maxwell element is the same, 
the model is best analyzed in terms of complex stiffness. Let the model con- 
sist of n Maxwell elements connected in parallel with a spring of stiffness k^, 
and a dashpot with damping coefficient c^. Also, let; the Jth element have a 
spring of stiffness kj and a dashpot with a damping coefficient of c^ . For 
solid materials, k^ will be finite and vrill correspond to the static stiffness 
and c^ may be zero. A relaxation time * measure of the time 

required for stress relaxation in the jth Maxwell element. A straightforward 
algebraic manipulation will show that the complex stiffness of the Maxwell 
iBodel is described by the function [9] 
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For brevity. Equation (Dl) is written as 
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If Kl and K2 are determined experimentally as a function of frequency, then 
Equations (D3) may be used to compute the various model coefficients. A least* 
squares method may be used to accomplish this task. Let the experimental data 
be tabulated at frequencies c. , 1 £ v ^ ">> where ro is the total number of data 

V 

points and also for generality let the functions Kl and K2 be weighted bv con- 
stants W1 and W2, then the sum of squared deviation (SSD) is defined as 
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where 
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Equation (D4) may be differentiated with respect to all the model coefficients 
and derivatives set equal to zero for minimum SSD. These resulting eqixatlons 
are to be solved simultaneously for the required coefficients. It is clear 
that the coefficients appear in a line-'r fashion but Tj or are nonlinear 
and are, therefore, determined by Iteratl. method. 

For any set of values for i s j s n, the linear coefficients k^, c^, and 

kj, 1 s j s n are determined by solving the ''normal equations" for a conventional 
linear regression analysis [l4]. 

Differentiation of Equation (D4) with respect to the nonlinear coefficient gives 
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For minimum SSD, the derivatives given by the above equations must be zero. 
Thus, the n simultaneous equations for 1 ^ j ^ n become 
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The above equations are solved by the Newton-Raphson Iteration method. The 
successive solutions are obtained by inverting the Jacobian matrix. 
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(D7) 


The elements of the Jacobian matrix are derived by straightforward differentia- 
tion of Equation (D6) 
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and h. 


The various derivatives of functions f and h are derivable by the relations 
defined in Equations 
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To summarize the iterative sdlution procedure the linear parameters, k^, c^, 
kj, 1 £ j £ n are solved for any set of iterates Oj > 1 ^ j ^ n. The new set 
of iterates Oj > 1 ^ j £ n are obtained by solving the linear algebraic Equa> 
tions (D7). 


It should be noted that such a least squares analysis may sometimes result in 
negative coefficients, which are physically meaningless. Under these conditions, 
it might seem reasonable to guess the possible relaxation frequencies from the 
trend of the data and determine the linear parameters by carrying out a 
simple linear regression. 

In order to Judge the reasonable values of Qj, 1 ^ J ^ n, the qualitative 

trends of Equation (d 4) roust be examined. It is dear that K1 is formed 

by linear superposition of f^ and K2 is obtained by superposing the functions 

hj. The functions f and h are plotted as functions of o)/Q in Figure d3. 

Thus, superposing f^ over a static stiffness k^ will result in an increasing 

K1 as a function of frequency. In the case of K2, some maxima or minima may 

be observed, depending on the values of Qj. If n * 1, then clearly K2 will 

peak at cj a If n > 2 and 0^^ and Q 2 quite far apart, then a minimum 

will be observed in K2 at cu somewhere between 0^^ and O 2 . If n a 2 and 0^^ and 

0^ quite close, Chen a maximum will be observed in K2 at cu somewhere be* 

tween 0^ and 02 * A non*zero value for coefficient c^ clearly indicates a linear 

increase in K2 with frequency. If no such trend ^s present in the data, c^ 

may be set equal to zero. These qualitative trends will establish the possible 

values of Q., 1 £ j s; n, and also might help determine the number of elements, 
j 

n, required in a model to fit the experimental data. 
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Fig. D3 Characteristic Functions Defining the Behavior of 
Visco-Elastic Models 




2. Voigt Model 

As shown tn Figure 1)2 let a generalized Voigt model be defined by n Voigt 
elements connected in series with a spring of stiffness k^. Also, let the 
jth element have a spring of stiffness and a dashpot with damping coeffi- 
cient Cj. If this model is applied to solid materials, and if the static 
spring with stiffness k^ is not present then the springs in all elements must 
have non-zero stiffness. This model is best analyzed in terms of complex com- 
pliance H(o)) which may be expressed as 
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(b + 
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1/k , b. = 1/k. and t. = c./k. 

O J j j j J 


(D12) 


The time in this case, is generally called retardation time since it is a 
measure of time required for the extension of the spring, in the Voigt element, 
to its equilibrium length while retarded by the dashpot. 


Equation (D12) for brevity may be expressed as 


H(o)) = H1 (cd) - i H2((o) 


where 
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HI (to) “ b + ^ b g (co) 
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and hj is same as defined in Equation (D4) , 
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It is clearly noted that Equation (01^)18 similar to Equation (D3) with = 0, 
In fact, If kjji kj. fj« 1 s j s n, and Kl^, K2^, 1 i v S ro in Equation (D3) 
are respectively replaced by b^, b^, g^, 1 s j s n, and Hl^, ^2^, 1 S v s m, 
th^n the resulting equation will be Identical to Equation (DlA). Thus, all 
thu analysis outlined for the Maxwell model Is also applicable in the case of 
the Voigt model. 


The derivatives of the function gj are easily obtained 
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All the qualitative features described in the case of the Maxwell model are also 
relevant here when the experimental data is plotted in terms of complex compli> 
ance. However, the function g^ must be used instead of f^. This is also plotted 
in Figure d3. 
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NOMENCLATURE 


2 

A Cross sectional area of elastomer specimen (m ) 

Input, shaker table amplitude (n) 

^2 Output, resonant mass amplitude (m) 

a Hc’^ht of elastomer specimen (ro) 

b Thickness of shear specimen (m) 

c Damping coefficient (N>sec/m) 

c^ Critical damping (see Equation Cl, Appendix C) (N-sec/m) 

c^ Elastomer damping coefficient (N-sec/m) 

Lower cylinder damping coefficients (N-sec/m) 

c^ Damping as defined in viscoelastic models in Appendix D (N-sec/m) 

c^ Upper cylinder damping coefficient (N-sec/m) 

Cj^,C 2 Variable damping coefficients as defined in Table 1 (N-sec/m) 

d Diameter of compression specimen (m) 

2 

E Young's Modulus of Elasticity (N/m ) 

F Static Force (N) 

f,g,h Characteristic functions defined in Appendix D 

2 

G Shear modulus of the material (N/m ) 

H Complex compliance (■ Hi — tH2) (m/N) 

i Imaginary unit ( ■ yTi) 

K Complex stiffness ( » Kl + iK2) (N/m) 

k Dynamic stiffness (N/m) 

k^ Elastomer dynamic stiffness (N/m) 

k^ Stiffness of lower cylinder (N/m) 

k^ Stiffness of mechanical spring (N/m) 


1 

I 

I 

I 

I 

i 

1 

I 

I 
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■f 

j 

i 

if 
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k Static stiffness defined in viscoelastic models in Appendix D 

(N/m) 

Stiffness of upper cylinder (N/m) 

k^ Stiffness of guide bearings (N/m) 

k ^ Static stiffness of elastomer specimen (N/m) 

St 

m Resonant mass (kg) 

m^cc Equivalent mass of accelerometer (kg) 

Equivalent mass of mechanical spring (kg) 

m . Equivalent mass of lower piston (kg) 

pi 

ffipu Equivalent mass of upper piston (kg) 

m Mass of the hub in the guide bearing (kg) 

8 

5 Shape factor 

8 Shore A hardness 

t Time (sec) 

u Dimensionless stiffness 

V Dimensionless damping 

W1,W2 Weights used in least squares analysis 
w Widtii of shear specimen (m) 

z Height of shear specimen (m) 

a Transmissibility ( ■ IA 2 /A 2 I ) 

3 Defined in Equation (C5) of Appendix C 

6 Static deflection (m) 

e Strain ( ■ 1 - X) 

Q Damping ratio as defined in Equation (Cl) of Appendix C 

\ Ratio of strained to unstrained length 

2 

a Stress (N/m ) 
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Relaxation or retardation time (aec) 

Phase angle (degree) 

Dimensionless frequency 

Hatural frequency of mechanical spring (rad/sec) 
Natural frequency (rad /sec) 

Natural frequency of guide bearings (rad/sec) 
Frequency (rad/sec) 
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